Higgsino dark matter model consistent with galactic 
cosmic ray data and possibility of discovery at LHC-7 



O 
(N 



Oh' 
I ■ 

Oh' 



> 
m 
o 



X 



Ning Chen/ Daniel Feldman,^ Zuowci Liu/ Pran Nath,"^ and Gregory Pcim'^ 

^C.N. Yang Institute for Theoretical Physics, Stony Brook University, Stony Brook, New York 11794, USA 
^Michigan Center for Theoretical Physics, University of Michigan, Ann Arbor, Michigan 48109, USA 
Department of Physics, Northeastern University, Boston, Massachusetts 02115, USA 

A solution to the PAMELA positron excess with Higgsino dark matter within extended super- 
gravity grand unified (SUGRA) models is proposed. The models are compliant with the photon 
constraints recently set by Fermi-LAT and produce positron as well as antiproton fluxes consis- 
tent with the PAMELA experiment. The SUGRA models considered have an extended hidden 
sector with extra degrees of freedom which allow for a satisfaction of relic density consistent with 
WMAP. The Higgsino models are also consistent with the CDMS-H and XENONIOO data and are 
discoverable at LHC-7 with 1 fb^^ of luminosity. The models are testable on several fronts. 

PACS numbers: 95.35.-l-d, 12.60.Jv, 13.85.Rm, 98.70.Sa 



INTRODUCTION 



Recently, experiments detecting galactic cosmic rays 
have begun to probe the nature of the dark matter in the 
halo. The large excess observed of high energy positrons 
in the PAMELA experiment [ij and the null results in 
the search for gamma ray lines with the Fermi-LAT ex- 
periment present a challenge for particle theory. Some 
particle physics explanations have been proposed to ex- 
plain the PAMELA data consistent with the relic abun- 
dance of dark matter including: a Breit-Wigner enhance- 
ment 0, a nonperturbative Sommerfeld enhancement 
0,01, and other possibilities (il-flB. A nonthermal cos- 
mological history is also a solution |lll4l3| . Several astro- 
physics explanations have also been sought [31 . Within 
supersymmetry the positron excess can arise from the 
annihilation of neutralinos [lightest supersymmetric par- 
ticle (LSP)l into W+W- and/or ZZ. This comes about 
when the LSP is a pure wino [ill [l3j , a mixed wino-bino 
[gI. [15| or a Higgsino [Tij - fioj . However, a wino LSP pro- 
duces a large amount of monochromatic photons in its 
annihilation products which is edging close to the cur- 
rent upper limit set by the Fermi-LAT data [3l ■ 

Here we present a supersymmetric model, which in 
contrast to other proposed models, has mostly a Higgsino 
LSP and can explain the relic abundance of dark matter. 
In addition, we show that such a model fits the positron 
excess from PAMELA [l| and is consistent with the an- 
tiproton flux, as well as with data from monochromatic 
photons that arise via loop diagrams in the neutralino 
annihilation processes xx ~^ 77,7^ [IS- We note that 
a bino-likc LSP can also explain the PAMELA positron 
data when a substantial size boost factor from the halo 
is allowed [2l|. 

The monochromatic photon constraints from Fermi be- 
come very relevant when one tries to fit the PAMELA 
positron data via dark matter annihilations in the galac- 
tic halo as the cross section needed to explain such data is 



much larger than the naive estimation of the dark matter 
annihilation cross section from a thermal history. This 
expectation for the relic abundance however can be mod- 
ified which will be discussed. Such a modification can 
open new parameter space in SUSY models where the 
relic density of dark matter is consistent with observa- 
tions and the flux of cosmics from dark matter at present 
temperatures can account for the data. This has im- 
plications for signatures of supersymmetry at the Large 
Hadron Collider in the frameworks we discuss below. 



EXTENDED ABELIAN MODELS AND 
ENHANCEMENT OF RELIC ABUNDANCE 

The simplest extension of the standard model (SM) 
which is gauge invariant, renormalizable, and unitary 

in. A 



arises through a Stueckclbcrg mechanism [22 



U{1) gauge boson V,,. gains mass M through a Stueck- 
clbcrg mechanism [2j] by directly absorbing an axion 
field a through the combination (MVJ^ -I- d^aY which 
is gauge invariant under the transformation SV^j, = dp,X, 
5a ~ —MX, and thus a transition to the unitary gauge 
produces a massive vector gauge boson without the ne- 
cessity of a Higgs mechanism. It is also well known 
that the Stueckelberg mechanism arises quite naturally 
from a Green-Schwarz mechanism [ist with appropriate 
transformations. Further, in reduction of higher dimen- 
sional theories the masses of the Kaluza-Klein states arise 
from a Stueckelberg mechanism and not from a Higgs 
mechanism. The Stueckclbcrg mechanism is indeed quite 
generic in string theories (see e.g. [1^), in extended su- 
pergravity theories, and in the compactification of higher 
dimensional theories (for a review see [27j). 

Interesting new physics arises if there is a hidden sector 
with minimally a U{1) gauge field that mixes with the 
hypercharge of the SM sector. A supersymmetric gener- 
alization of the Stueckelberg mechanism leads to an ex- 
tended neutralino sector, i.e., where for each extra U{l)x 



factor one has two extra Majorana fields (Stinos) which 
mix with the minimal supersymmetric standard model 
(MSSM) neutralinos. The above considerations general- 
ize to a set of Abehan U{l)x gauge groups, and such 
extensions lead to a mixing between fields in each sector 
via gauge kinetic energy mixings and mass mixings. 

We implement this extension to study a class of su- 
pergravity unified models which allow the possibility of 
explaining the PAMELA data without recourse to large 
clump factors in the halo of the Galaxy. We uncover a 
new situation where the LSP is actually a nearly pure 
Higgsino under radiative electroweak symmetry break- 
ing with mass in the range ~ (110 — 190) GeV with the 
hidden sector components of the LSP being suppressed. 

Thus we consider a supergravity grand unified model 
[28I [2^ having an extra hidden sector with a product 
gauge group U{l)x which mixes with the hypercharge 
via mass terms generated by the Stueckelberg mecha- 
nism and without loss of generality via gauge kinetic 
mixing. For simplicity, we give a summary for the case 
of a single U{l)x, and the generalization for a prod- 
uct gauge group follows analogously. In the vector sec- 
tor the mass mixing and gauge kinetic energy mixing 
is of the form -2MjfMyX''y^ - {S/2)X^"'Yf,^ and in 
the ncutralino sector the mass mixing is of the form 
'4'stiMxXx+MYXY) + h.c. while the kinetic mixing leads 
to —iS{Xxcr ■ dXy + <-> X)), where X denotes the hid- 
den sector U{1) and Y is the hypercharge of the MSSM; 
ipst is a fermonic field that arises out of a chiral Stueck- 
elberg supermultiplet and My ■ My and 6 are small, 
i.e. on the order of 10~^ or smaller [23|. Such additional 
states remain in contact with the thermal bath prior to 
freezeout in the early universe. In the absence of hidden 
sector soft masses, a direct study of the mass matrix in 
the ncutralino sector gives rise to a mass degeneracy g^d 
for the hidden sector neutralinos with the LSP, which in 
turn, can have a degeneracy (7vis with other visible sec- 
tor sparticles [sot . Coannihilations can then produce an 
enhancement of the relic density by a factor /e @ so 
that 
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Generalizing to the case of a J7(l)^ extended hidden sec- 
tor (7hid = 2n, and thus for the case ^vis = 1, one finds 
Je = (2n + If which gives Je = 25(49) for n = 2(3). 
In this extended model the ncutralino mass matrix will 
be (4 -I- 2n) x (4 -I- 2n) dimensional. We assume that the 
LSP lies in the visible (MSSM) sector. Because of coan- 
nihilations in the visible sector, the full enhancement is 
never achieved, however one finds large enhancements of 
size (10 — 20) or larger with a degenerate hidden sec- 
tor and only 2 — 3 additional t^(l)s which is sufficient 
for compatibility with the WMAP constraint [sH since 
the models considered have the relic density in the range 
~ (2 — 6) X lO"'^ if there were no Abelian hidden sector. 



Model 


mo 


Ml 


M2 


M3 


^0 


tan /3 


1 

M 


1 

Ml 


M2 


M3 


PI 


1033 


1600 


1051 


120 


2058 


13 


195 


683 


836 


259 


P2 


1150 


1600 


1080 


160 


2080 


15 


152 


684 


859 


347 


P3 


950 


1425 


1820 


748 


1925 


25 


109 


617 


1453 


1589 


WB 


2000 


400 


210 


200 


300 


5 


562 


170 


163 


441 



TABLE I: Parameters which produce an LSP which are 
mostly Higgsino (P1-P3), or mixed wino-bino, WB. Here 
"lo(^o) is the universal scalar mass (trilinear coupling), 
Ml , M2 , M3 are the gaugino masses at the GUT scale for the 
gauge groups U{1)y, SU{2)l, SU{3)c and tan/3 is the ratio of 
the two Higgs vacuum expectation values in the MSSM. The 
parameters that enter the neutralino mass matrix at scale 
Q = i/M~M^ are (/i', Mi, M2, M3), where fi' is the Higgs 
mixing parameter. The models have also been run through 
both SuSpect and SOFTSUSY via micrOMEGAs [H. Here 
mP„t=173.1 GeV. 





Einasto 


NFW 


Isothermal 


Model 


/ \ theory 


180[190] 4.4[2.3] 
130[150] 5.3[2.5] 
90[110] 4.3[0.7] 
150[160] 5.9[2.0] 


6.1[3.2] 
7.3[3.5] 
6.0[1.0] 
8.2[2.7] 


10.4[5.5] 
12.6[6.0] 
10.3[1.7] 
14.1[4.7] 


PI 
P2 
P3 
WB 


0.24[0.08] 
0.23[0.09] 
0.18[0.09] 
7.00[1.29] 



TABLE II: Cross sections {av)-,z and {av}~f-f upper lim- 
its (10~^^cm^/s) 0] for 3 halo profiles (Einasto, Navarro- 
Frenk- White (NFW), and Isothermal) along with predictions 
for (P1-P3) and WB. The mostly Higgsino models (P1-P3) 
are unconstrained by any profile while the mixed wino-bino 
model, WB, is on the edge. 



LOW MASS HIGGSINO LSP IN EXTENDED 
SUGRA AND FERMI PHOTONS 

We discuss now the details of the Higgsino-like neu- 
tralino models. Since the extra weak mixing discussed 
above is small, it has negligible effects on the soft pa- 
rameters at the weak scale. The model parameters that 
dictate annihilation cross sections can then be described 
by the input parameters given in Table (|l]). The mod- 
els (P1-P3) listed in Table(jT| have a neutralino that is 
dominantly a Higgsino, with about 2% remaining in the 
gaugino content. For comparison, we also exhibit the 
mixed wino-bino (WB) model @ which has a significant 
wino content ^ 49% of the total eigencontent along with 
a comparable bino content. All four models satisfy the 
current experimental constraints from flavor physics and 
limits on sparticle masses (see e.g. [s^). Their neutralino 
masses lie in the range (110-190) GeV and have a spin 
independent cross section of size (5-10) x 10~^^cm^ con- 
sistent with the upper bounds from the CDMS-II and 
XENONIOO [131. Further, some of the models possess 
several rather light sparticles in their spectra, namely 
the charginos, neutralinos, gluino and in some cases the 
stop, and are thus good candidates for discovery at the 
LHC. 

In Table (|IH ) we give the theoretical predictions of the 
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Higgsino LSP models for the and 77 modes and 
exhibit the current upper limits from the Fermi-LAT 
search for photon lines using three different halo pro- 
files. One finds that the theoretical predictions for the 
Higgsino models (P1-P3), are well below the current up- 
per bounds from Fermi-LAT, by about a factor of 10, for 
the most restrictive profile, while the mixed wino-bino 
model, WB, is close to the edge of the limits. There 
are sources of photons arising from bremsstrahlung that 
could mimic the line signature of monochromatic pho- 
tons. The contributions from bremsstrahlung to the line 
signals can be significant or even dominant over the ones 



from the loop processes [35|]. However, the additive ef 



fects from bremsstrahlung to the line source are small 
for the models considered here which have the dark mat- 
ter in the mass range ^ (110 — 190) GeV. This is 
due to the fact that the maximal energy the photon 
can carry is E^'^'^^ = M^{1 - M^^/M^) in the process 
XX ~^ WW J, and the energy of the monochromatic pho- 
tons via XX iX is = M^[l - M'^/{AM^)]. Thus 
the photons arising from the process XX ~^ WW^ have 
energy whose location in the energy spectrum is at least 
^ (23 — 40) GeV below the monochromatic photons in 
7Z, 77 final states for dark matter mass in the range 
~ (110 — 190) GeV. For a related discussion see [36| . 

The Higgsino models typically have a small \i and large 
Too and lie on the boundary of the radiative electroweak 
symmetry breaking curve, i.e., the Hyperbolic Branch 
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38| . It is the smallness of [i relative to the soft gaug- 
ino masses that makes the three lightest particles, the two 
lightest neutralinos and the lighter chargino, essentially 
degenerate in mass (ssj . In this region /i (and some of 
the sparticle spectrum) is very sensitive to small changes 
in the input parameters at the GUT scale. On the other 
hand since [i is small, one is in a less fine tuned region. 
Alternately, instead of working down from the high scale, 
one could simply generate these Higgsino-like LSPs di- 
rectly by inputs at the weak scale. We have checked this 
for the models discussed here. This is evident from Ta- 
blejU. 



POSITRONS FROM HIGGSINOS 
AND MIXED WINOS 

Next, we discuss the positron excess prediction in the 
Higgsino-like model. In Higgsino and wino models, the 
high energy positron flux can arise from WW and ZZ 
production from the neutralino annihilation in the halo 
with approximate cross sections at leading order [lij 
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where V = (W,Z), xy = M^/M^, Cw = 16 (1) for 
Higgsino (wino) models and the ZZ production is only 



significant for Higgsino models where Cz = 32 cos^(6'iy). 
For the models (P1-P3) the LSP is mostly a Higgsino 
with only a very small portion being gaugino. Here the 
cross sections that enter in the positron excess are size 
{av)ixX^WW,ZZ)u 



iggsino 



< 4 X lO-'^'^cm'^/s 



The positron flux from the Higgsino dark matter can 
be described semianalytically (for early work, see (isj). 
The flux enters as a solution to the diffusion loss equation, 
which is solved in a region with a cylindrical boundary. 
The particle physics depends on (crw)haio, and dN/dE, 
the fragmentation functions / energy distributions ISl. 
The astrophysics depends on the dark matter profile 39| , 
and on the energy loss in the flux from the presence of 
magnetic fields and from scattering off galactic photons. 
A boost factor which parametrizes the possible local in- 
homogeneities of the dark matter distribution can be 
present. Recent results from N-body simulations indi- 
cate that large dark matter clumps within the halo are 
unlikely [i^ [4l[. The boost B we consider here is small, 
as low as ~ (2 — 3). The background taken is consistent 
with the GALPROP [13] model generated in Ref. 1 of 
13|. The antiproton flux follows rather analogously (for 
an overview and some fits see e.g. (isj). In this analysis 
the antiproton backgrounds are consistent with [3] , and 
the results for the pure wino case considered are consis- 
tent with 

The full analysis is exhibited in the upper left panel 
of Fig.(IT) where we show fits to the PAMELA positron 
fraction Ij . For comparison we also show the essentially 
pure wino case, which requires no boost (clump), but as 
mentioned in the introduction, will generally lead to an 
overproduction of photons. Model P3 requires a boost of 
only - (2 - 3) as the LSP is hght, - 110 GeV. For this 
case, the p flux is slightly larger at lower kinetic energy, 
but still consistent with the data. A pure wino at 110 
GeV would give a cross section about 10 times larger rel- 
ative to the Higgsino model at 110 GeV. Including the 
boost factor of 3 for the Higgsino model, the pure wino 
is then (3 — 4) times stronger in its flux, and this is an- 
other reason a pure wino at 110 GeV would fail - it would 
overproduce the antiprotons, whereas the Higgsino with 
minimal boost is consistent. Thus, in the upper right 
panel of Fig.([T]) we give a comparison of the^ flux with 
the recently released data from Ref. 3 of [l|. Indeed 
it is seen that the theoretical prediction of the p flux is 
in perfectly good accord with this data. We note there 
are other processes beyond the leading order that could 
produce SM gauge boson final states. For diboson final 
states, these corrections are rather small and lead to a 
small shift downward in the clump factor used (see Ref. 
2 of [4^). The minimal boost utilized here is rather dif- 
ferent compared to those in analyses of bino-like LSPs 
which use boosts of size 10^ or larger 2l| to fit the data. 
The analysis we present does not attempt to explain the 
high energy e + e data H, 47 1 . This could be explained 
with an additional electron sourcell3l. 
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SUGRA models 

PAMELA Data 
PAMELA Data (Stat) 




^— Higgsino P3, B = 


3,M = 


108 GeV 


Higgsino P1,B = 


5,M = 


190 GeV 


— I — Higgsino P2, B = 


5,M = 


148 GeV 


Pure Wino, B = 


1,M = 


180 GeV 


mixed WB, B = 


2,M = 


162 GeV 


Background Fit 







PAMELA Data 2010 
Signals + Background 



E GeV 




Kinetic Energy T GeV 



^ 80r 



- 7 ToV □ Background 
■■IS - I lev p^^^ ^3^3 g^yj 

P2 (m,=450 GeV) 

WB (m,=591 GeV) 



p(j)>150GeV - 




400 



600 800 1000 1200 1400 1600 1800 

Effective Mass (GeV) 



SM Background 
PI (m =348 GeV) 

P2 (m =450 GeV) 

g 

WB (m,=591 GeV) 




200 



p^(j^) > 1 50 GeV 
Pj%,y,)^'iO GeV 
^^>^50 GeV 
S,^ > 0.2 
n(l) = 
e {e, n) 



400 600 800 1000 1200 

di-Jet invariant Mass (GeV) 



FIG. 1: (Color online) Top left: PAMELA positron excess and the Higgsino models (P1-P3). The wino dominated model is 
also shown for comparison along with a mixed wino-bino model, WB. Top right: The PAMELA p flux and the predictions are 
seen to be compatible with the data. Equal dark matter densities and boosts are taken in both top panels. Lower left: SUSY 
plus background events vs effective mass at 1 fb~^ for the signature cuts shown in the panel for models PI, P2 and WB. Lower 
right: SUSY plus background events vs the di-jet invariant mass at 1 fb"'^ for signature cuts shown in the panel for the models 
PI, P2 and WB. Both lower plots do not show model P3 due to its suppressed cross section at LHC-7. The legends labeling 
the model names (not colors) indicate the model class in all figures. 



SIGNATURE ANALYSIS: LHC, ^ = 7 TEV 



As mentioned above, some of the colored sparticles in 
the Higgsino-like models are rather light which is encour- 
aging for possible early discovery of this class of models 
at LHC — 7 [4^. To achieve a significance necessary for 
discovery, i.e., S > maxjsV^, 10}, it is essential to have 
a reliable SM background computation. In our anal- 
ysis we simulate the SM backgrounds [i^ using Mad- 
Graph 4.4 [H^l for parton level processes, PYTHIA 6.4 
for hadronization and PGS-4 for detector simulation (sij . 
The ^-tagging efficiency in PGS-4 is based on the tech- 
nical design reports of CMS and ATLAS ^ (see 
The sparticle spectrum and branching ratios for the sig- 



nal analysis is generated u sing computational packages 
for supersymmetric models (32|. 

The models we consider for the LHC-7 analysis have 
rather light gluinos in the mass range (^ 350 — 600) GeV 
(see also [53[). The production cross sections for these 
models are dominated by gluino production and the 
branching fractions are dominated by either the radia- 
tive decay of the gluino g — gxl, g — >■ gx2 (Higgsino-like 
model PI) or a combination of the radiative decays above 
and the three body decays g — >■ xti^t + ^^-C-) (Higgsino- 
like model P2) or effectively just the 3 body decays pro- 
ducing both xt and X2 with substantial rates (mixed 
wino-bino model WB) . The subsequent decays follow 
from the chargino and neutralino into standard model 
quarks and leptons. Decays into the degenerate hidden 
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sector particles near the LSP mass are suppressed. 

In the lower left panel of Fig.(IT]) we give an analy- 
sis for the Higgsino models PI and P2 with the number 
of SUSY events in 40 GeV bins at 1 fb^^ of integrated 
luminosity vs the effective mass defined to be the sum 
of the pt of the four hardest jets plus missing energy. 
The cuts used are exhibited in the panel. In the lower 
right panel of Fig.([T]) we give an analysis for the Higgsino 
models PI and P2 with the number of SUSY events in 
50 GeV bins at 1 fb""'^ of integrated luminosity vs the 
di-jet invariant mass where the cuts used are exhibited 
in the panel. For comparison we also give an analysis of 
the mixed wino-bino model, WB, in both lower left and 
lower right panels. Since the gluino is relatively light 
and the squarks are heavier, the 3 body decays of the 
gluino dominate resulting in rich di-jet signals and effec- 
tive mass. We note that while the model P3 provides 
a good fit to the PAMELA data and its photon flux is 
an order of magnitude below the current limits, it has 
a heavy (~ 1.5 TeV) gluino and would not produce an 
identifiable signal in the early LHC data. 

CONCLUSION 

We have presented here a solution to the PAMELA 
data and the Fermi photon data with a Higgsino-likc 
LSP which can also be made compliant with WMAP. It is 
shown that the models considered are consistent with the 
current very stringent limits on 77 and production 
from Fermi-LAT which put the pure wino LSP models 
close to the edge of the upper limit of experiment. Fur- 
ther, the Higgsino LSP models are consistent with the 
upper limit from the XENONIOO experiment and will be 
testable in improved dark matter experiments. We find 
that LHC-7 can realistically probe these models up to 
gluino masses of ^ 600 GeV with 1 fb~^ of data. How- 
ever, one would need larger integrated luminosity to carry 
out precise mass reconstructions. The above presents 
an interesting possibility of having a low mass gluino 
from the radiative breaking of the electroweak symmetry 
which can be produced at the LHC in early runs and also 
having a mostly Higgsino LSP giving rise to PAMELA 
positron excess. Thus the class of models discussed here 
can be tested on multiple fronts. 
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